Decoherence as attenuation of mesoscopic echoes in a spin-chain channel 
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An initial local excitation in a confined quantum system evolves exploring the whole system, 
returning to the initial position as a mesoscopic echo at the Heisenberg time. We consider a two 
weakly coupled spin chains, a spin ladder, where one is a quantum channel while the other represents 
an environment. We quantify decoherence in the quantum channel through the attenuation of the 
mesoscopic echoes. We evaluate decoherence rates for different ratios between sources of amplitude 
fluctuation and dephasing in the inter-chain interaction Hamiltonian. The many-body dynamics is 
seen as a one-body evolution with a decoherence rate given by the Fermi golden rule. 
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I. INTRODUCTION 

Control of quantum dynamics is essential to achieve 
quantum information processing. Examples of this are 
implementations of quantum algorithms [1, 2] and quan- 
tum communications [3], |4| . The system involved in each 
of these processes interacts with an environment which 
degrades the quantunr correlations [5] . This loss of infor- 
mation, called decoherence, represents the main obstacle 
to achieve an efficient quantum processing. Thus, in or- 
der to avoid decoherence, it is mandatory to understand 
its processes. Many approaches exist to study decoher- 
ence in systems composed by few qubits [6-9]. However, 
real computer implementations that involve large num- 
ber of qubits require the dev elop ment of new approaches 
to characterize decoherence |10l4l2| . 

Recently, man y e fforts were done to characterize the 
quantum noise [13l4l5l | that produces decoherence in 
quantum channels. These channels connect two quantum 
systenrs enabling the information transfer between them. 
Spin chains can be used to achieve this goal for short dis- 
tance comnrunications [16] avoiding interfaces between 
the static system and the information carriers. Whereas 
all these studies were oriented to pure-state communi- 
cation processes, experimental realizations of pure-state 
dynamics are a major challenge [17]. Alternatively, im- 
plementations of quantunr computation with NMR can 
not deal with pure states but have to resort to statistical 
mixtures of spin ensembles instead [18-20]. State transfer 
in spin ensembles was done in a ring of spins with many- 
body interactions in the solid state [2l|, [22j • An initial 
local polarization propagates around the ring. The return 
of the initial excitation is evidenced through the construc- 
tive interference that reappears at the Heisenberg tinre 
tjj ~ h/A, with A the typical mean energy level spacing, 
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as a form of polarization revival called the mesoscopic 
echo [23]. There, the polarization amplitude and phase 
of a given nuclear spin within the ring was monitored as 
a function of time. However, the many body nature of 
spin-spin interactions strongly compromises an optinrum 
transfer. Thus, a nruch more efficient transfer was ob- 
served through the experimental implementation of an 
effective XY Hamiltonian (i.e. flip-flop processes) in a 
spin chain [24] by exploiting the J-coupling in the liq- 
uid phase, where dipolar interaction becomes negligible. 
The main reason for this result is that the many-body 
dynamics of an XY Hamiltonian is mappable to a one- 
body evolution [25]. Thus, in this case, the Heisenberg 
time is proportional to the system size M instead of the 
2 M value that shows in a complex many-body Hamilto- 
nian. In Ref. [2J], the evolution of the initial excitation 
was monitored in all the spins of the quantum channel. 
Comparisons with theoretical calculations showed the ef- 
fect of decoherence manifested in the attenuation of the 
interference intensities which is nrore evident in the de- 
cay of the mesoscopic echoes. This experimental break- 
through enabled various theoretical proposals for perfect 
state transfer |26l . \2l\ . More recently, implementations 
of spin-chains were done in solid-state NMR [28| by im- 
plementing a Double Quantum Hamiltonian (i.e. flip- 
flip+flop-flop processes) which, in turn, is mappable to 
an XY interaction [29]. However, the progress toward 
the proposal to observe the mesoscopic echoes of these 
systems [301 was experimentally limited by length inho- 
mogeneity [3l|. 

In this work we propose to use the attenuation of the 
mesoscopic echoes in a spin chain as a sensor of the de- 
coherence produced by an uncontrolled spin bath, in this 
case a second spin chain. The spin-spin interaction within 
each chain is given by an XY Hamiltonian where the en- 
semble dynamics can be solved analytically |32h34| ] . Once 
the chains are laterally coupled to form a spin ladder, the 
quantum dynamics becomes truly many-body and the 
analytical solution is no longer possible. Moreover, nu- 
merical solutions are difficult to obtain as a consequence 
of the exponential increase of the Hilbert space dimen- 



sion [35|. We show that, within certain range of the ra- 
tio between the inter-chain and intra-chain interactions, 
the evolution of a local excitation in the many-body sys- 
tem (spin ladder) can be obtained as a one-body dynam- 
ics (isolated chain evolution) plus a decoherence process 
given by a Fermi golden rule (FGR) . We characterize the 
decoherence rate for different kinds of inter-chain inter- 
actions by controlling the contribution of the different 
sources: pure dephasing (Ising interaction) and ampli- 
tude fluctuations (XY interaction). Thus, we show that 
this auto-correlation function becomes an effective and 
practical tool for studying decoherence. Moreover, the 
exploration nature of the mesoscopic echo provides the 
sender (Alice) with the information of the decoherence 
effects in the quantum channel by contrasting with the 
expected result of the isolated dynamics. This could en- 
able the implementation of a convenient error correction 
code without needing of a receiver. 

In the next section, we introduce the system and 
present the numerical calculations of the local polariza- 
tion evolution for different inter-chain interactions. From 
these results we obtain the decoherence rates from the 
attenuation of the mesoscopic echoes. In the third sec- 
tion, for a solvable small system, we characterize analyti- 
cally the decoherence rates. Instead of using the quantum 
master equation which is the most standard framework 
ado pted to describe the system-environment interaction 
[36l |37||. we use the Keldysh non-equilibrium formalism 
which leads to an integral solution of the Schrodinger 
equation [33, 34, 38, 39]. Finally, we discuss the conclu- 



II. SPIN DYNAMICS IN A NON-ISOLATED 
SPIN-CHAIN CHANNEL 

A. System 

Let us consider the polarization dynamics in a spin 
system composed by two parallel spin chains transversely 
coupled. The topology of the interaction network forms 
a ladder. Chain I represents the quantum channel com- 
posed by M spins, and chain II is another M spin system 
that perturbs the quantum transfer in chain I, see Fig. 
[IJa). The spin Haiuiltonian is given by 



'Htotai — ~Hi + Hn + Ht, 



(1) 
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represents the cv-th spin chain Hamiltonian that takes 
into account the XY interaction between neighbor 1/2- 
spins within the chain. The spin chains interact through 



the following transversal Hamiltonian 
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where by controlling the ratio a/b one determines the 
nature of the interaction. The polarization evolution is 
considered for different transversal interactions: XY in- 
teraction (a/b = 0); isotropic (Heisenberg) interaction 
(a/b = 1); and truncated dipolar interaction (a/b = —2), 
which are typical in NMR experiments [3a, [37|. Be- 
sides these three typical NMR interactions we will also 
consider the evolution under two additional transversal 
Hamiltonian with the following parameters Hi (a/b = 2) 
and H<2 (a/b = 3.6). We will focus on regular systems 
where all longitudinal couplings within chains are taken 

equals J n +\ „ = Jx ( n =/" M). The same is applied to the 
transversal couplings J„ = J y . We choose these param- 
eter conditions because the interference effects are more 
pronounced and the mesoscopic echo degradation can be 
easily evaluated. 



a) 



b) 
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Figure 1: (Color online): (a) Schematic representation of the 
spin ladder system given by the Hamiltonian of Eq. JTJ . Each 
spin chain has M = 5 spins, (b) A system like in (a) but 
with periodic boundary conditions between spins 1 and M. 
Light colors represent the quantum channel. Marked spin 
(red online) represents the site where the initial condition is 
placed. 



B. Decoherence characterization based on the 
attenuation of mesoscopic echoes 

In order to study the mesoscopic echoes of a spin en- 
semble, we calculate the evolution of a local polarization 



within chain I (quantum channel) through the spin auto- 
correlation function 133, 34], 



Pi.At) = 



(y cq \sr(o)srm\y 



(G) 



cq / 



This gives the local polarization in the z direction on site 
1 at time t providing that the system was in its equi- 
librium state plus a local excitation on site 1 at time 



t = 0. Here, S[ a) (t) = e 



iH t , 



t/ft£(") z e -m t „tait/fi is the 



spin operator in the Heisenberg representation and \^ e q/ 
is the many-body state corresponding to thermal equi- 
librium, that is a mixture with appropriate statistical 
weights, of all possible states with different number N of 
spins up. In the regime of NMR spin dynamics k^,T is 
much higher than any energy scale of the system. Then, 
all the statistical weights can be taken equal [3a, l37| . 
As numerical method we alternate between the standard 
diagonalization by spin projection subspaces [21] or the 
Quantum Parallelism algorithm [35|, which involves the 
evolution of a few superposition states representing the 
whole ensemble. Although the last is much more efficient 
for larger samples, for the considered system sizes, both 
involve similar computing time. As the results coincide, 
no further comment is devoted to this point. 

By assuming J y equal to zero, the local polarization 
dynamics at spin 1 in chain I is equivalent to the lo- 
cal evolution in an isolated spin chain, P 1 s ° latcd (£). This 

function is defined with Eq. © replacing "Htotai by "Hi 



Since the XY spin- 



and S[ I)Z (t) = ^Hjt/h S m' e -w^/\ 
spin interaction within the chain only couples nearest 
neighbors, the system dynamics in the high tempera- 
ture regime can be obtained analytically as a one-body 
dynamics [33, [3J|. Upper line (black) in Figure [2]^a) 
shows the solution of pj s °i ated (£) f or a S pi n chain with 
M = 5. We can observe the presence of the fc-th meso- 
sco pic echo at a time t^ E proportional to the chain size 
M [2l|. As J y increases, the degrees of freedom of chain 
II start modifying the observed dynamics of pj s ° lated (£) 
and, in general, no analytical solution is known. To 
obtain the polarization dynamics Pi j i(i) of this many- 
body problem, we solve numerically the time dependent 
Schrodingcr equation by using the Trotter-Suzuki decom- 
position |40j. Color curves of Fig. HJa) show the local po- 
larization Pi,i(t) for different values of J y corresponding 
to an isotropic transversal Hamiltonian. They evidence 
the degradation of mesoscopic echoes proportionally to 
J y . Fig. \T[i>) shows Pi,i{t) for different forms of the 
transversal interaction Hamiltonian for a fixed value of 
J y . Due to the XY term in the transversal Hamiltonian, 
the polarization is transferred back and forth between 
chain I and II. This process, plus the dephasing induced 
by the Ising contribution, produce the progressive spread- 
ing of the polarization among all the spins thus degrading 
the strong recurrences. This is manifested in the mean 
value of the local polarization at longer times, t 3> H/ J x , 
that approximately tends to 1/ (2M). Thus, one observes 
a decrease of the mesoscopic echoes with evolution time 



as well as a gradual increase of the background polar- 
ization at times between echoes. We do not consider an 
Ising transversal interaction (b/a = 0) since it does not 
involve transfer of polarization between chains leading to 
a different value of the mean local magnetization at long 
times. In the latter case, the mean local polarization at 
long times tends to a value that is close to \jM instead of 
the 1/ (2M ) value. This difference in the final state avoid 
a direct comparison with the polarization dynamics from 
Hamiltonians that contain an XY term requiring extra 
manipulations. 




10 20 30 40 50 60 70 80 90 

time [ft/Jx] 



Figure 2: (Color online): Polarization at site 1, Pi,i(t), of the 
spin ladder of Fig[TJa). (a) Upper (black) curve represents 
the dynamics of an isolated chain given by p^° lated (£) . Color 
curves correspond to different values of the coupling J y for a 
transversal isotropic interaction between chains, (b) Compar- 
ison of the dynamics behavior between different transversal 
interaction natures (XY, Isotropic, dipolar) by keeping fixed 
the coupling ratio J y /J x = 0.1. 

In order to characterize decoherence of the quantum 
channel, we measure the attenuation of the mesoscopic 
echoes as compared with the local polarization in the iso- 
lated spin-chain channel, P 1 s ° latod (£). Figure[5]shows that 
the ratio P 1A {tf E )/Pi s f at * d (tf E ) decreases as a func- 
tion of t^f E . Within the regime |J y /J x | < 1 we ob- 
serve that this ratio has an exponential dependence of the 
form exp[— t^ E /t^]. This exponential decay of the initial 
polarization cannot hold for all times until P l 1 (t) =0 
since the system size is finite and the total polariza- 



tion is conserved within the system. Thus, for a later 
time the correlations within chain II are manifested at 
site 1, and the ratio P hl (tf E )/Pff ated (tf E ) evidences 
a complex behavior instead of an exponential one. The 
numerical results show that the crossover between the 
two temporal behavior appears at times proportional to 
£r ck t^ ln[ J^t^/Zi] [4l|. To characterize the exponen- 
tial decay of the mesoscopic echoes we calculate the ratio 
• p i.i( i r E )/- p i?i lated (^ E ) for several values of tf E previ- 
ous to this crossover time. The accuracy of this charac- 
terization improves proportionally to the order k of the 
mesoscopic echo. Thus, to avoid dealing with too long 
times or large number of spins, we increase the order of 
observable echoes without changing the physical behav- 
ior of the system by putting periodic boundary condition 
in Wij\. This "closes" the spin chain into a ring as shown 
in figure QJb) . This choice approximately duplicate the 
number of mesoscopic echoes for a fixed evolution time 
with respect to the observed in the chain. 

Since chain II does not have infinite degrees of free- 
dom it does not represent a reservoir. Thus, the presence 
of recurrences in the quotient P h i(tf E )/ P^^itf®) is 
expected instead of an exponential law. In figure [3J the 
function -log(P 11 (i)/Pi s f lated (i)) is plotted for different 
values of J y as a function of time. Due to the dynamics 
within the channel each curve oscillates around zero. The 
maximum values are associated to the temporal region 
where pj s ° lated (f) > P 1A (t), that is, the region where the 
mesoscopic echoes are manifested. The minimum values, 
on the other hand, correspond to the temporal regions 
between mesoscopic echoes, where p^ s ° latod (£) < P 11 (i). 
It can be seen that an exponential law is identified in 
the envelop of the peaks that appears for each curve (see 
dashed line) . In figure [4] are plotted the values of the 
maxima of figure [3] for J y = 0.09J X and J y = 0.10J X as 
a function of time. The peaks start to get apart from its 
envelop at different times proportional to £r depending 
on the intensities of the transverse interaction. This be- 
havior reflects the recurrences due to the finite nature of 
chain II, which becomes relevant for earlier times as the 
intensities J y increases. In order to quantify the expo- 
nential behavior corresponding to each value of the trans- 
verse interaction J y , we fit the maximum values shown 
in Fig. [3] to a linear function. Thus, we extract the deco- 
herence time r^ from the slope of the curve for different 
values of the transversal interaction J y . The number of 
points (number of echoes) used in this procedure varies 
depending on the value of J y as can be seen in Fig. 21 
Besides, the point (0,0) is not included in the fit since 
for very short times t <C h/J x the behavior of quantum 
dynamics is not exponential but it starts quadratic with 
t. 

The existence of an exponential law in the decay of 
the mesoscopic echoes where r^ ex \HJ X /J 2 \ plus the fact 
that its time regime is bounded by a time proportional to 
tn ex T0 In [J x T,p/h] suggests that the decay of mesoscopic 
echoes in chain I could be described by the Self Consistent 
Fermi Golden Rule (SC-FGR) analyzed in Ref. [U. 
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Figure 3: (Color online): Natural logarithm of the ratio 
Pii(t)/Pi s f latod (t) for different values of the transversal cou- 
pling J y as a function of time. The interaction between chains 
has an isotropic nature. The dashed line is a guide to the eye. 
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XY (a/b = 0) 


(0.96 ± 0.04) %f 


Isotropic (a/b — 1) 


(1.47 ± 0.05) %£ 


dipolar (a/b — —2) 


(2-7 ±0.1) %L 


Pi (a/b = 2) 


(0.81 ± 0.04) %f 


P 2 (a/b = 3.6) 


(1-8 ± 0.1) J^ 



Table I: Proportional constants between l/t^> and \J y \ 2 /(hj x 
for different interactions between chains. 



Figure O shows the obtained values of r</, as a function 



of 



''yl 



evidencing a linear dependence. The curves cor- 
respond to transversal interaction of different nature by 
varying the ratio a/b in Eq. J3]): XY (a/b — 0), isotropic 
or Heisenberg (a/b = I), truncated dipolar (a/b = —2), 
Hi (a/b = 2) and H 2 (a/b = 3.6). The slopes of the 
linear curves are shown in table [I] 

The linear behavior observed in figure [5] confirms that, 
even when we are dealing with closed systems in which 
the energetic spectrum of the environment (chain II) is 
not continuous, the characteristic time r^ agrees with 
the one dictated by the FGR. The conditions needed to 
derive the FGR in molecular spin systems was already 
addressed in ref. |4lJ. There, the behavior of the numer- 
ical simulations on clusters with less than 20 spins for 
times shorter than £jj was observed to coincide perfectly 
with the theoretical results. 

In the present case, 5+5 spins, the unperturbed Hamil- 
tonian is represented by Pi +TLu, while the perturbation 
is Ht- The initial state is not an eigenstate of %i + Tin, 
but instead it is a superposition of the 2 5 eigenstates of 
chain I. Since the initial state is spatially localized, all the 



0.6 



0.4 



CL 0.2 
0.0 



J /J =0.09 

y x 



10 20 30 40 50 60 

time [h/J v ] 

^0.8] J /J =: „10 
y x 

I 0.6 

Q_ 
>^ 0.4 



0.2 



0.0 



10 20 30 40 50 60 

time [ fi I J ] 



Figure 4: (Color online): Maximum values of Fig. [3] for two 
different ratios J y /J x as a function of time. The slope of the 
linear fit (solid line) gives the value of r^ associated to the 
perturbation J y . 



eigenstates of chain I contribute with equal weight to this 
superposition. Therefore, the decay of a local polariza- 
tion involves a sort of average over all possible initial and 
final states. This justifies the use of the expression for 
the characteristic time given by the Fermi Golden Rule 

El 
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2lT 



H T 



N , 



(7) 



where 



H T 



is a characteristic value of the coupling in- 
teractions between chains that has to be determined and 
No represents the density of directly connected states, i.e. 
states of the states of chain II connected with the initial 
state through the perturbation Ht- This interaction con- 
sists of two processes: the flip-flop (or XY) interaction 
and the Ising interaction weighted by the factors a and b 
respectively. Thus, the evolution of the initial localized 
state generated in chain I decays due to the coupling with 
the "environment" represented by chain II. Equation ([7]) 
involves the assumption that No is similar for the XY 
and Ising interaction, i.e., iVg CY = N™ = Nq. This den- 
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Figure 5: (Color online): Decoherence rates — for different 

transversal Hamiltonians as a function of \J y \ /hj x . The lin- 
eal dependence indicates the good agreement with the FGR 
behavior. Different slopes indicate that the decoherence rate 
is composed by two contributions: One is due to the injection 
of polarization produced by the XY term of Ht and the other 
is associated to the dephasing process given by the Ising term. 



sity of states can be approximated as the inverse of the 
second moment of Tin, and then No oc 1/J X . The cross 

_ 2 , , 

terms proportional to ab are canceled out in "Ht [36| , 

thus, under these approximations each interaction terms 
between chains becomes independent of each other ob- 
taining 
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r XY 



B*W 



fiJx 



(8) 



(9) 



where A and B represent constants associated with each 
interaction. We determine these constants by using the 
numerical results of table |TJ Equation ((SJ can be rewrit- 
ten as 



yxY 



A 
B : 



J ZZ 



1 
B 



(10) 



where x zz = a 2 T<j,\J y \ / (HJ K ) and y XY = 
b 2 T<f,\J y \ J (HJ X ). Figure ^ shows the curve y XY 
vs. x zz (squares) for the three standard NMR Hamil- 
tonians, XY (b = 1, a = 0), Isotropic (b = 1, a = 1), 
dipolar (b — I, a = —2) and for other two arbitrary 
Hamiltonians Tii and Tii with parameters (b = 0.5, 
a = 1) and (b — 0.5, a = 1.8) respectively. By doing a 
linear fitting (solid line), we obtain 



A = (0.49 ± 0.08) and B = (1.00 ± 0.06) 



(11) 
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III. TWO-SPIN CHANNEL COUPLED TO A 
SPIN BATH: ANALYTICAL SOLUTION 

A. System 

In this section, we consider the simplest quantum spin 
channel: a two spin system that could act as a SWAP 
gate. Each of these spins is coupled with an independent 
spin environment as is sketched in Fig. [7] Since the sys- 
tem describes a Rabi oscillation that could be damped by 
the interaction of the spin bath, it contains all the essen- 
tial ingredients of the spin channel case, with the further 
advantage that it could be solved analytically and to asses 
in detail how the environment disturbs it. Therefore, it 
allows the exact characterization of the decoherence rate 
in order to compare with the numerical results of the 
previous section. 

The Hamiltonian of the two spin system is given by Hi 
of Eq. © with M = 2. The system-environment Hamil- 
tonian is equivalent to Eq. ^ with M = 2. Finally, 
the environment is represented by two independent and 
semi-infinite linear chains whose Hamiltonians are given 
by Tin of Eq. © with M — >■ 00. By neglecting the inter- 
action between the semi-infinite portions we ensure that 
no correlation could appear. We call them Hel and Her, 
see figure [7J As the theoretical framework we adopt the 
Keldysh formalism |4a| in a form that was developed for 
electronic excitations [44, 45] and then extensively devel- 
oped to solve the spin dynamics in pre sence of system- 
environment interactions |33l . [3J, |38|, (39fl ■ More recently, 
this formalism also proved to be useful to address thermal 
transport |46j |. 

In this work, we briefly discuss the main points of the 
formalism required to arrive to the solution. We start by 
establishing the relation between spin and fermion oper- 
ators at a given site n by applying the Jordan- Wigner 
transformation (JWT) [25( 



-0.5 0.0 0.5 1.0 1.5 2.0 



s n = c^ exp{i7r 



(12) 



Xz 



Figure 6: (Color online): j/ XY vs. x zz for different natures 
of the transversal interaction, where x zz — a 2 r^ \J y \ / (HJ X ) 
and 2/xy = b 2 i> |J y | 2 /ftj x 
fitting to the points. 



The solid line is given by a linear 



Here, c^ and c n stand for the creation and destruction 
fermionic operators, and S^ are the rising and lowering 
spin operator, S„ = S* ± iS^. Within the fermionic de- 
scription, the spin-up and spin-down states correspond 
to an occupied and not-occupied fermionic state respec- 
tively. 



It is important to note that the role of the interaction 
nature on the decoherence tinre is manifested through the 
constants A and B. The origin of this factors will be clar- 
ified in the next section by using the Keldysh form of the 
quantum field many-body theory. There, a microscopic 
model is applied to a similar spin system. We will obtain 
analytical expressions that agree well with the numerical 
values obtained here. 
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Figure 7: (Color online): Schematic representation of the two 
spin channel coupled to independent spin environments. 
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After applying the JWT to the total two spin Hamil- 
tonian 



^2-spin = fi l+ fa 



H 



EL 



Wer, 



(13) 



it is possible to rewrite the different contributions of 

#2-spin ^ ig given by 

%i = V 12 (c+c 2 + c+c,) , (14) 

where the hopping amplitude between states 1 and 2 is 

Vm = J/2 

and the transition between them occurs at the natural 
Rabi frequency 



cj = 2V 12 /h. 
The environment is represented by 



(15) 



n 



EL 



%ER 



E^0 + -* 

i=l 



L, i'-h,i+l 



C L,i+l C L,i 



oo 



C R,i+l + C R,i+l C R, 



(16) 



electron in the first site of the left (right) reservoir cor- 
respond to an Ising interaction between spins 



V zz 



1L 

/ZZ 
^R 



V ZZ 

> / 91 



aJlL : 

a^2R 



aJ y 
: a J v 



Analogously, the XY (flip-flop) component of the inter- 
chain interaction is associated to the hopping amplitudes 






1/ 2R 



bJy/2. 



(19) 



Note that the Ising term of the system-environment inter- 
action in the spin problem is not completely analogous to 
the Coulomb interaction within the fermionic description, 
since after the JWT the exchange term of the Coulomb 
interaction is not present. The last three terms in each 
of the brackets of Eq. (fT8|) do not involve an interaction 
between the system and the environment. Thus they just 
modify the potential energies of sites 1, 2 and site 1 of 
both environments. Since the potential energy can be 
controlled externally to maximize the polarization trans- 
fer in a NMR experiment [47I . |48| . these terms will be 
neglected for the present calculations. 



B. Spin dynamics within the Keldysh formalism 



where cj ^(cl j) is the creation (destruction) fermionic op- 
erator for the environment that interacts with the system 
state 1 and similarly for Cr j(c Ri ) which belongs to the 
environment that interacts with the system state 2. Vl 
and Vr stand for the hopping term between neighboring 
sites within each of the environments. They correspond 
to the XY (flip-flop) interaction along the x direction. 



y L = V R = J x /2 



(17) 



Finally, the transversal or system-environment interac- 
tion takes the form 



% 



V = ^1L (CIA + C+C L ,l) 

.^ 1. 



V, 



ZZ 
1L 



■-.+ 



- 



C L,1 C L,1 C 1 c l — 2 C L,1 C L,1 ~~ 2 Cl ° l 4 



l/ 2R 



V ZZ 



where 



1 



C R,1 C R,1 C 2 C 2 2 C R-' iCr > 1 



l.+ „ 1 

2 2 2 4 



(18) 



V& Y = bJ 1L /2 



V, 



XY 
2R 



6J 2 r/2 



are the hopping amplitudes, due to the XY term of the 
interaction, between the site 1 and 2 in the quantum 
channel with the left and right environments respec- 
tively. The standard direct integral of the Coulomb in- 
teraction of an electron (fermion) in state 1 (2) with an 



Even at room temperature k^T is much higher than 
any energy scale involve in an NMR experiment. The lo- 
cal polarization given by Eq. ^ can be expressed within 
the Keldysh formalism as [33] 



*m,ny') 



i "m,mr! V 



1. 



(20) 



Here, G^ m (t, t) is a particular case of the general particle 
density function [49| 



G< (*a,«i) = i (*n. c .|^(ii)c n (t 2 ) |*„.e.) 



(21) 



The initial polarized state is described by the non- 
equilibrium state l^n.e.) — cj I'J'cq.) formed by creating 
an excitation at t = on the n-th site. Therefore, the 
non-equilibrium density G^^t^t) of Eq. (|20|l depends 
implicitly on the index n that indicates the site of the ini- 
tial excitation. The expression for this initial condition 
can be expressed as 



G< ; (0,0) = ^ (6 k ,i + 2AP S k , n 5 n ,i) , 



(22) 



where the first term describes the equilibrium density 
which is identical for all sites and does not contribute 
to the dynamics. The second term represents the non- 
equilibrium contribution where only for the n-th site is 
different from zero. Notice that this initial state repre- 
sents a non-correlated initial state, where the only ele- 
ments different from zero are those satisfying k — I. The 
factor AP accounts for the excess of excitation at site 
n and is the responsible of the observed dynamics. Its 
value ranges from zero (lack of excitation) to a maximum 
of 1/2. 
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The initial density function, Eq. (|22)) . evolves under 
the Schrodinger equation that could be expressed in the 
Danielewicz form [43| given by the expression 



G< (t 2 , h) = h z G n (t 2 , 0) G< (0, 0) G A (0, ii) 



di fc d^G R (t 2 ,t k ) £< (t k , ti) G A (t u ti) . (23) 



Here, G < (t%, ti) is a 2 x 2 density matrix whose elements 
G^ n (t 2 ,ti) are restricted to TO, n € {1,2}. Similarly, 
G R (i2,^i) represents an effective evolution operator in 
this reduced space whose elements, the retarded Green's 
function 

G* n (t 2 ,*i)=K m (t 1 ,t a )] t (24) 

= 9(t 2 ,ti) [G mn (t 2 ,h) - G mn (t 2 ,ti)], 

(25) 

describe the probability of finding an excitation at site 
to after it was placed at site n and evolved under the 
total Hamiltonian for a time t 2 — t\ . The injection self- 
energy, 11 < , takes into account the effects of the envi- 



ronment. The first term of Eq. (j23j) stands for the "co- 
herent" evolution since it preserves the memory of the 
initial excitation, while the second term contains "inco- 
herent reinjections" described by the injection self-energy 
that compensates any leak from the coherent evolution 
d. 

In absence of Ht, the retarded Green's function for 
the system is easily evaluated in its energy representation 
given by the following expression 



G UH (e) = / G UK (t) exp[iei/ft]di = [el - Hi] 



(26) 



Conversely, with presence of T-Lt, the interacting Green's 
function defines the reduced effective Hamiltonian 



H cff .(e) = £l-[G R ( £ )] ^Hi 



E*(e) 



(27) 



and the self-energies £ R (e) [5fJ|, where the exact per- 
turbed dynamics is contained in the nonlinear depen- 
dence of the self-energies E R on e. For infinite reservoirs, 
ReS R (e zy ) = Re (e u — e°) represents the "shift" of the 
system's eigen-energies s° with e v the eigen- energies of 
the interacting problem. The imaginary part of E R , 



2 Im S R (e„) /h = l/r T = 2T T /h, 



(28) 



accounts for their "decay rate" into collective system- 
environment eigenstates in agreement with the Self- 
Consistent Fermi Golden Rule [41], i.e. the evolution 
with H c ff. is non-unitary 

We use a perturbative expansion on Ht to build up ex- 
pressions for the particle (hole) self-energies 'S < ^ > '(tk, t{) 
as well as for the retarded (advanced) self-energies 
X R ( A )(£fc, ti). Under the wide band assumption (or fast 



fluctuation approximation), where the dynamics of exci- 
tations within the environments are faster than the rel- 
evant time scales of the system Vl, Vr 3> V12, we ob- 
tain for the decay rates of the XY and Ising system- 
environment interaction [34j, l38l . |39( | the following expres- 
sions 



nXY 



2tt 



— 1 IT — — — V-i 



IT 



iXY 



1L 



2tt 



l 2T _ \ V r 



XY 

2R 



2 1 
2 1 



nV R 



and 



^zz 



2tt 



fc 1 it - -T- \y 1 



ZZI 



API 



2 r zz 

n 



2T 



n 

2tt 

h 



3tt 2 U l 



— \v zz \ 



AP|) 



3tt 2 Ur 



(29) 
(30) 

(31) 
(32) 



Here, it is considered that the left and right environments 
are in equilibrium. Thus, the occupation for any of their 
sites is fL,R = (1/2 + APi^), with AP1.2 representing 
the occupation excess on the left and right environment 
respectively. Due to the wide band approximation the 
decay rates are time and energy independent, and in con- 
sequence H e ff. does not depend on e 



H e g. = 




XY 



where Tit = T 1T 



XY 



■ rf| and r 2T = r 2T 



rff . In order 



to obtain a comparison with the numerical results of the 
previous section, we assume that the excess of occupa- 
tion in the left (right) environment APi (AP2) is very 
small (APi, AP2 <C ^), condition that is well satisfied in 
the high temperature regime, and that the hopping am- 
plitudes satisfy V^ = U 2 Z R Z , V$? = Vg? and V h = Vr. 
They ensure that the decay rates to the left and right 
environments are identical, i.e. Tit = T2T = Tt- Un- 
der these conditions, the propagator has a simple de- 
pendence on t given by G R (t) = G 0R (t) e~ rTt ^ h , where 

G?f(*) = G%f(t) = H C0S (T*) and G ™(t) = G "?(i) = 
jrsm(tft). Then Eq. $F5§ becomes 



,-t/TT. 



G< (t, t) = h? G UH (t) G< (0, 0) G UA (-*) e 

f dt t G m (t - U) S< (U) G 0A (ti - t) e-C*-**)/'*, (33) 
Jo 

which is a generalized Landauer-Biittiker equation [44j, 
|45j. This equation is complemented with the injection 
self-energy S < , which takes into account for particles 
that return to the system after an interaction with the 



environment, and is expressed as |34 1381. 13' 

i (\ 



and 






2Tl z hG< + h2T^ Y - - + APi 



XY 



2T^hGf % + h2T% 



h\2 
i 1 



H2 



dij, 



h (34) 
(35) 



where in the last expression we used the assumption 
APi,AP 2 <l/2. 



C. Characterization of the decoherent processes 



We solve Eq. (f3"3"]l together with the injection self- 
energy of Eq. (|34p with an initial condition given by an 
excitation of AP on site 1, i.e. fG<- (0,0) = APS a Si :r 
For this purpose we follow the strategy used in Refs. [3a . 
HI]. Replacing Eq. ([Ml) into Eq. J33]) and identifying the 
interaction rate of Eq. (|28p , we get two coupled equations 
for Gfi and G^ 2 given by 



§G< U (t,i) = fi 2 

(22) 



G u ft (t) 



APe- t/TT 



G°S (t-ti) 

(21) 



G°?2 (t-ti) 

(22) 



e -(t-t 4 )/r T [A E < (,.)] di . 



-(t-t«)/Tr rziv< 



[f S< (U)} Mi. (36) 



The first term is the probability that a particle, initially 
on site 1, is found at time t on site 1 (or 2) having survived 
the interaction with the environment. The second and 
third terms describe particles whose last interaction with 
the environment, at time U, occurred at site 1 and 2 
respectively. Noticeably, in the first term of Eq. (|36|) the 
environment, though giving the exponential decay, does 
not affect the frequency of the two spin system given 
in G°S (t). Modification of to requires the dynamical 

(21) 

feedback contained in the next terms. The solution of Eq. 
(|3"6"]l involves a Laplace transform. This solution, for the 
general case where the system-environment interaction 
involves both XY and Ising terms, gives for the local 
polarization of Eq. ([2"0"]) the following expression 



Pn(t) = 2AP 
1 



l e -2T™t/H 



+ - 



2 cos ((/>) 
where 



[( W + h ? )£ + 0] e -( 2r * Y+r T Z ) i / r > 



^-(sy 



-\/(£) - 1 



hu > r| z 
hujo < rf z 
huj > r zz 
huj < rf z 



(37) 

(38) 
(39) 



pZZ 

1 T 



tan(^) = --J- 

huj 



(40) 



From these expressions it is possible to determine the 
observable frequency u) and the decoherence time as 
the slowest of the two competing interaction rates: 
ty [(21^ + 1^) ±7,]. 

It is interesting to remark that the effect of the lateral 
chains on the two spin system coupled through an Ising 
system-environment interaction can produce observables 
with non- linear dependences on Ht- We find a non- 
analyticity in these functions enabled by the infinite de- 
grees of freedom of the environment [5lJ (i.e. the ther- 
modynamic limit). Here, they are incorporated through 
the respective imaginary part of the self-energy, ?i/r zz , 
i.e. the FGR. Hence, the non-analyticity of uj and r^ 
on the control parameter fiwo/r zz at the critical value 
Swo/Tt 2 = lj indicates a switch between two dynamical 
regimes. In previous works, we identified this behavior as 
a Quantum Dynamical Phase Transition [38l . 1481 ] . This 
can be interpreted as a disruption of the environment 
into the dynamical nature of the system, a form of the 
Quantum Zeno Effect (QZE), which states that quantum 
dynamics is slowed down by a frequent measurement pro- 
cess (52j. Very recently, this dynamical phase transition 
was evidenced as a particular scale invariance in the fluc- 
tuations in the number of photons emitted by a driven 
two level system |53| . However, here we choose to work 
in the regime where decoherence rate is still weak to pro- 
duce such transitions. The evaluation of the conditions to 
observe such dynamical phase transition or a cross-over 
in the dynamics dimensionality J54| is an open problem 
that deserves further studies. 

For a pure Ising system-environment interaction, Tt = 
i.e., r^ Y = 0, we recover the expression found in 
Ref. [38]. For a pure XY system-environment inter- 
action Tt = r^ Y , i.e., r zz = 0, it is noticeably that 
the observable frequency and decoherence time depend 
linearly on %t- In particular, the oscillation frequency 
coincides with that of the isolated system, i.e. uj = uq. 
This is due to the symmetry of the decay rates to both 
environments, Tit = r 2 T- In absence of this symmetry 
this effect is not observed as can be seen for a sp in system 
interacting with only one environment [34, 39, 48]. 

It is remarkably that in the limit Hujq 3> T zz , the so- 
lution (J37]) tends to the solution of a two spin system 
only coupled to the environment through one spin of the 
system [55j. Thus, within this regime it is impossible to 
identify whether the two spin system is coupled to two 
wide band environments or to a single one. 

From Eq. (13"T|) it is straightforward to obtain the 
decoherence time within the regime J y <C J x where 



rzz 
i T 



l/i> = (r£ z + 2r£ Y ) jh. In order to compare with the 
numerical results of previous section we identify the intra- 
chain and inter-chain hopping as 



Vl = V R = J x /2 



(41) 
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and V{1 



XY 



V, 



XY 
2R 



bJy/2, 



(42) 



while the through space inter-chain Coulomb coupling 
results: 



V zz 



,/ 2R 



aJ v 



(43) 



Thus, taking in mind the high temperature regime 
(APi, AP2 < |) , we obtain for the decay rates r zz and 



r^ Y the following values 



-D 



XY 



and 



T^ZZ _ 

J- T — 



277 1 u , 



2tt 



2 1 
1 ^J x 



ft ' aJyl 3tt 2 J x 



(44) 



(45) 



It is interesting to note that these two decay rates are 
not equal for an isotropic system-environment interaction 
(a/b =1) [331 where its ratio is given by 



rzz o 

X T _ « . 849 _ 



r* Y 



3tt 



(46) 



Moreover, increasing the occupation of the environments 
APi and AP2 the ratio is reduced. This result con- 
trasts with the usual assumption of taking them equal 
[H, Szl, S [H, lint Thus, here we show that always 
an XY system-environment interaction is more effective 
than the Ising one to destroy coherence. Observing the 
decoherence rate 1/r^ = (ri| z + 2r* Y ) /H, an extra fac- 
tor of 1/2 reduce the Ising decoherence rate. In addi- 
tion, one can decrease the Ising decoherence process by 
increasing the occupation within the reservoir. 

By using Eqs. (l4"4l and (J45I) we obtain for the deco- 
herence rate 



1A* = 



1 



1 



_ZZ ' „XY 

I A. I A 



XY\ 



T^ + 2T^ Y )/h 



4 a 2 \J y 
3n hJ„ 



b 2 Wy 

UJy- 



0.424 



a 2 \J y \ 

hJ K 



b 2 \J y \ 



(47) 

(48) 
(49) 



These analytical results are in notable agreement with 
those obtained from the numerical solutions of the dy- 
namics of spin systems computed in the preceding sec- 
tion, Eqs. (JHJ) and (TJTJI . In particular it confirms 
a notable effect repeatedly observed experimentally: a 
stronger interaction along the chain results in a weaker 
decoherence [31|. Indeed, Eq. U7] states that a fast in- 
chain dynamics makes the already slow inter-chain dy- 
namics even slower. This is a form of the QZE that is 
manifested experimentally in the spin diffusion in low- 
dimensional crystals. Slightly different crystals showed 
an unexpected dimensional crossover as a function of a 
structural parameter [5j| . This crossover was described 
as a QZE where the internal degrees of freedom act as 



measurement apparatus [54j . The concept that the mea- 
surement is played by an interaction with another quan- 
tum object, or simply another degree of freedom of the 
subsystem investigated, was independently and fully for- 
malized by recasting it in terms of an adiabatic theorem 
by Pascazio and collaborators J53|. Conversely, a strong 
inter-chain interaction can even lead to a freeze of the 
SWAP dynamics as described by Eq. [35] and fully char- 
acterized experimentally |48l |. 

Getting more into more precise details, the Keldysh 
description of the spin dynamics enabled to obtain from 
first principles the origin of the coefficients A and B ob- 
tained in previous section on numerical grounds. While 
the precise scale depends on the success of the model to 
yield a representative local density of directly connected 
states, the relation between A and B depends on the na- 
ture of the XY and Ising interaction and hence has a 
universal meaning 



IV. CONCLUSION 

We have studied numerically the quantum evolution of 
a localized initial excitation within a spin-chain channel 
weakly coupled (\J y /J x \ < 1) to a second lateral spin- 
chain, i.e. a spin ladder system. We characterized the 
decoherence rate of the channel by measuring the atten- 
uation of mesoscopic echoes with respect to the isolated 
channel evolution for different kinds of the inter-chain in- 
teractions. We showed that the XY system-environment 
interaction is more effectively to destroy quantum coher- 
ences than the Ising system-environment interaction. No- 
tably, by increasing the environment occupation, APi^, 
one can reduce even more the Ising decoherence rate. 

The decoherence characterization was possible by re- 
sorting to the analytical solution of a two spin channel, 
where each of the spins are coupled to independent spin 
environments in a fast fluctuation regime. That was done 
by using the Keldysh formalism. We showed that a quan- 
tum dynamical phase transition only appear when an 
Ising system-environment interaction is present. Thus, 
for an XY system-environment interaction the bare two- 
spin oscillation frequency holds. The very good agree- 
ment between the decoherence rates of the spin-chain 
channel and the two spin system confirms that, within 
the weakly coupled regime, the finite lateral chain be- 
haves as an environment in a fast fluctuation regime. 
Within this regime, in consequence the complex many- 
body evolution of the spin ladder can be obtained from a 
one-body dynamics plus an exponential decoherence pro- 
cess dictated by the Fermi Golden Rule. Although sim- 
ple in nature, this statement has important experimental 
consequences. It helps to explain in terms of the QZE a 
notable experimental observation in low dimensional spin 
dynamics: the stronger interaction along the chain dom- 
inates over the weaker inter-chain inducing a stability of 
one-dimensional dynamics against perturbations by spins 
outside the chain [3l], [58] . A natural suggestion is that in 
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a quantum channel one can use the attenuation of meso- 
scopic echoes as a tool for characterizing decoherence in 
the channel without using a receiver at the other end. 



E.P.D. thank the Alexander von Humboldt Foundation 
for a Research Scientist Fellowship. P.R.L. and H.M.P. 
are members of the Research Career of CONICET. 



Acknowledgments 

We acknowledge support from Fundacion Antorchas, 
CONICET, FoNCyT, and SeCyT-UNC. G.A.A. and 



[1] I. L. Chuang, N. Gershenfeld, and M. Kubinec, Phys. [28 

Rev. Lett. 80, 3408 (1998). 
[2] L. M. K. Vandersypen et al., Nature 414, 883 (2001). [29 

[3] M. Riebe et al., Nature 429, 734 (2004). 
[4] M. D. Barrett et al., Nature 429, 737 (2004). [30 

[5] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 
[6] F. H. L. Koppens et al., Nature (London) 442, 766 [31 

(2006). 
[7] Y. Wu and X. Li and L. M. Duan and D. G. Steel and 

D. Gammon, Phys. Rev. Lett. 96, 087402 (2006). [32 

[8] M. Grajcar et al., Phys. Rev. Lett. 96, 047006 (2006). 
[9] M. Riebe et al., Phys. Rev. Lett. 97, 220407 (2006). [33 

[10] H. G. Krojanski and D. Suter, Phys. Rev. Lett. 93, 

090501 (2004). [34 

[11] H. G. Krojanski and D. Suter, Phys. Rev. Lett. 97, 

150503 (2006). [35 

[12] G. A. Alvarez and D. Suter, to be published in Phys. 

Rev. Lett. (2010). larXiv: 1004.50031 [36 

[13] Michael A. Nielsen, Isaac L. Chuang, Quantum Compu- 
tation and Quantum Information, Cambridge University [37 

Press (2000). 
[14] D. Burgarth and S. Bose, Phys. Rev. A 73, 062321 

(2006). [38 

[15] J. M. Cai, Z. W. Zhou, and G. C. Guo, Phys. Rev. A 74, 

022328 (2006). [39 

[16] S. Bose, Phys. Rev. Lett. 91, 207901 (2003). 
[17] A Quantum Information Science and Technology [40 

Roadmap, http://qist.lanl.gov/ (2004). 
[18] D. G. Cory, A. F. Fahmy, and T. F. Havel, Proc. Natl. 

Acad. Sci. USA 94, 1634 (1997). 
[19] N. A. Gershenfeld and I. L. Chuang, Science 275, 350 [41 

(1997). 
[20] E. Knill and R. Laflamme, Phys. Rev. Lett. 81, 5672 [42 

(1998). [43 

[21] H. M. Pastawski, P. R. Levstein and G. Usaj, Phys. Rev. [44 

Lett. 75, 4310 (1995). [45 

[22] H. M. Pastawski, G. Usaj, and P. R. Levstein, Chem. [46 

Phys. Lett. 261, 329 (1997). 
[23] V. N. Prigodin, B. L. Altshuler, K. B. Efetov and S. Iida, [47 

Phys. Rev. Lett., 72, 546 (1994). 
[24] Z. L. Madi, B. Brutscher, T. Schulte-Herbruggen, R. [ 

Briischweiler and R. R. Ernst, Chem. Phys. Lett. 268, 

300 (1997). [49 

[25] E. H. Lieb, T. Schultz, and D. C. Mattis, Ann. Phys. 

N.Y. 16, 407 (1961). [50 

[26] M. Christandl, N. Datta, A. Ekert and A. J. Landahl, 

Phys. Rev. Lett. 92, 187902 (2004). [51 

[27] P. Karbach and J. Stolze, Phys. Rev. A 72, 030301 

(2005). [52 



P. Cappellaro, C. Ramanathan, and D. G. Cory, Phys. 

Rev. A 76, 032317 (2007). 

S. I. Doronin and E. B. Fel'dman, Solid State Nucl. Magn. 

Reson. 28, 111 (2005). 

P. Cappellaro, C. Ramanathan, and D. G. Cory, Phys. 

Rev. Lett. 99, 250506 (2007). 

E. Rufeil-Fiori, C. M. Sanchez, F. Y. Oliva, H. M. 

Pastawski, and P. R. Levstein, Phys. Rev. A 79, 032324 

(2009). 

E. B. Fel'dman, R. Briischweiler and R. R. Ernst, Chem. 

Phys. Lett. 294, 297 (1998). 

E. P. Danieli, H. M. Pastawski, and P. R. Levstein, Chem. 

Phys. Lett. 384, 306 (2004). 

E. P. Danieli, H. M. Pastawski, and G. A. Alvarez, Chem. 

Phys. Lett. 402, 88 (2005). 

G. A. Alvarez, E. P. Danieli, P. R. Levstein and H. M. 

Pastawski, Phys. Rev. Lett. 101, 120503 (2008). 

A. Abragam, The principles of nuclear magnetism, 

Clarendon Press, Oxford (1961). 

RR. Ernst, G. Bodenhausen, A. Wokaun, Principles of 

nuclear magnetic resonance in one and two dimensions, 

Oxford University Press, Oxford (1987). 

E. P. Danieli, G. A. Alvarez, P. R. Levstein and H. M. 

Pastawski, Solid State Commun. 141, 422 (2007). 

G. A. Alvarez, E. P. Danieli, P. R. Levstein and H. M. 

Pastawski, Phys. Rev. A 75, 062116 (2007). 

H. D. Raedt and K. Michielsen, |quant-ph/04 06210 

Handbook of Theoretical and Computational Nanotech- 

nology. Quantum and Molecular Computing, Quantum 

Simulations (American Scientic Publishers, 2006). 

E. Rufeil Fiori and H. M. Pastawski, Chem. Phys. Lett. 

420, 35 (2006). 

P. Facchi and S. Pascazio, Physica A 271, 133 (1999). 

P. Danielewicz, Ann. Phys. 152, 239 (1984). 

H. M. Pastawski, Phys. Rev. B 44, 6329 fl991). 

H. M. Pastawski, Phys. Rev. B 46, 4053 (1992). 

L. Arrachea, G.S. Lozano, and A. A. Aligia, Phys. Rev. 

B 80, 014425 (2009). 

S. R. Hartmann and E. L. Hahn, Phys. Rev. 128, 2042 

(1962). 

G. A. Alvarez, E. P. Danieli, P. R. Levstein and H. M. 

Pastawski, J. Chem Phys. 124, 194507 (2006). 

L. V. Keldysh, Sov. Phys. JETP 20, 1018 (1965)[Zh. 

Eksp. Theor. Fiz. 47, 1515 (1964)]. 

P. R. Levstein, H. M. Pastawski and J. L. D'Amato, J. 

Phys.: Condens. Matter 2, 1781 (1990). 

S. Sachdev, Quantum Phase Transitions, (Cambridge U. 

P., 2001). 

P.B. Misra and E. C. G. Sudarshan, J. Math. Phys. 18, 



12 



756 (1977). 
[53] J. P. Garraham and I. Lesanovsky, Phys. Rev. Lett. 104, 

160601 (2010). 
[54] H. M. Pastawski and G. Usaj, Phys. Rev. B 57, 5017 

(1998). 
[55] G. A. Alvarez, P. R. Levstein and H. M. Pastawski, Phys- 

ica B 398, 438 (2007). 
[56] L. Miiller, A. Kumar, T. Baumann and R. R. Ernst, 

Phys. Rev. Lett. 32, 1402 (1974). 



[57] A. K. Chattah, G. A. Alvarez, P. R. Levstein, F. M. 

Cuchietti, H. M. Pastawski, J. Raya, and J. Hirschinger, 

J. Chem. Phys. 119, 7943 (2003). 
[58] P. R. Levstein, H. M. Pastawski, and R. Calvo, J. Phys.: 

Condens. Matter 3, 1877 (1991). 
[59] P. Facchi and S. Pascazio, Phys. Rev. Lett. 89, 080401 

(2002), and references therein. 



